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Abstract. Big Bang Nucleosynthesis is the theory of the production of 
the light element isotopes of D, 3 He, 4 He, and 7 Li. After a brief review of 
the essential elements of the standard Big Bang model at a temperature of 
about 1 MeV, the theoretical input and predictions of BBN are discussed. 
The theory is tested by the observational determinations of the light ele- 
ment abundances and the current status of these observations is reviewed. 
Concordance of standard model and the related observations is found over 
a limited range of the baryon-to-photon ratio, the single true parameter 
of the standard model. Implications of BBN on chemical evolution, dark 
matter and constraints on particle properties will be also discussed. 

1. Introduction 

Big Bang Nucleosynthesis (BBN) is the theory explaining the origins of the 
light elements D, 3 He, 4 He, and 7 Li and their primordial abundances. The 
theoretical framework for BBN is quite simple. It is based on Friedmann- 
Lemaitre-Robertson- Walker cosmology [1] and a network of nuclear reac- 
tions. We can further specify that the standard BBN model refers to ho- 
mogeneous nucleosynthesis in the context of a FLRW Universe with an 
electroweak standard model particle content, which for the purposes of 
BBN really amounts to assuming the existence of three nearly massless 
and nearly stable neutrinos. The predictions of BBN for the abundances of 
the light elements are determined by running a code which incorporates the 
nuclear network in a thermal (and cooling due to the expansion of the Uni- 
verse) bath. These predictions are then be compared with the observational 
determinations of the abundances. 
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In contrast to the theoretical side of BBN, the status of the observa- 
tional data has changed significantly in the last several years. There is more 
data on 4 He and 7 Li, and data on D and 3 He that was simply non-existent 
several years ago. For the most part, the inferred abundances of 4 He and 
7 Li have remained relatively fixed, giving us a higher degree of confidence 
in the assumed primordial abundances of these isotopes as is reflected in 
their observational uncertainties. Indeed, the abundances of 4 He and Li 
alone are sufficient to probe and test the theory and determine the single 
remaining parameter in the standard model [2], namely, the baryon-to- 
photon ratio, rj. In contrast, D and 3 He are highly dependent on models of 
chemical evolution ( 3 He is in addition dependent on the uncertain stellar 
yields of this isotope). New data from quasar absorption systems, on what 
may be primordial D/H is at this time disconcordant, different measure- 
ments give different abundances. As a consequence of the uncertainties in 
D and 3 He, one can hope to use the predictions based on 4 He and 7 Li in 
order to construct models of galactic chemical evolution. These results also 
have important implications for the amount of (non)-baryonic dark matter 
in the galaxy and on the number of allowed relativistic degrees of freedom 
at the time of BBN, commonly parameterized as N u . 

1.1. STANDARD MODEL BASICS 

Since one of the main inputs of the theoretical side of BBN is the standard 
hot big bang model, it will be useful to review some of the key concepts 
as they pertain to BBN. The metric for the FLRW model of the Universe 
is specified by two quantities, the curvature constant k, and the expansion 
scale factor R(t). At early times, the curvature is unimportant as can be 
seen from the Friedmann equation for the Hubble parameter 



where A is the cosmological constant. Since the density p scales as either 
R~ 3 (for a matter dominated universe) or R~ 4 (for a radiation dominated 
universe) , this term dominates over either the curvature or the cosmological 
constant. I will ignore both in what follows. 

The critical energy density p c is defined such that p = p c for k = 




(1) 



Pc = 3H 2 /8ttG n 



(2) 



In terms of the present value of the Hubble parameter this is, 



Pc = 1.88 x l(T 29 /i 2 gcm" 3 



(3) 
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where 

h Q = H / (100kmMpc~ V 1 ) (4) 
The cosmological density parameter is then defined by 

n = (5) 

Pc 

in terms of which the Priedmann equation, Eq. (1), can be rewritten as 
(with A = 0) 

(Q - 1)H 2 = A ( 6) 

so that k = 0, +1, — 1 corresponds to Q = 1, $7 > 1 and < 1. (Very) broad 
observational limits on h and Q are 

0.4 < h a < 1.0 0.1 < n < 2 (7) 

The value of Q, at least on relatively small scales, seems to depend on 
scale. Indeed, the contribution to Q from visible matter associated with 
stars and hot gas is quite small, Q ~ 0.003 — 0.01. On somewhat larger 
scales, that of galactic halos or small groups of galaxies, f2 ps 0.02 — 0.1. 
On galaxy cluster scales, it appears that Q may be as large as 0.3. And 
while there is some evidence, the observations are far from conclusive in 
indicating a value of £1 as large as 1. It is however possible to obtain a 
bound on the product, Q.h 2 from 



H t = / {i-n + n/x)- 1/2 dx 

Jo 



(8) 



(for A = 0). For t > 12Gyr, and Q < 1, fM 2 < 0.3 (This is true even if 
A/ O.i 

As indicated above, BBN takes place during the radiation dominated 
epoch which lasts roughly to the period of recombination (somewhat earlier 
when dark matter is included) which occurs when electrons and protons 
form neutral hydrogen through e~ + p — > H +7 at a temperature Tr~ few 
xlO 3 K ~1 eV. For T < Tr, photons are decoupled while for T > Tr, 
photons are in thermal equilibrium. Today, the content of the microwave 
background consists of photons with T Q = 2.728 ± .002 K [3]. The energy 
density of photons in the background can be calculated from 



Pi 



J E y dn y (9) 



where the density of states is given by 



dn y = -%[exp{E y /T) - l]' 1 q 2 dq (10) 
2ir z 
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and g 7 = 2 is the number of spin polarizations for the photon, E 1 = q is just 
the photon energy (momentum). (I am using units such that Ti = c = ks = 
1 and will do so through the remainder of these lectures.) Integrating (9) 
gives 

P, = yT* (11) 

which is the familiar blackbody result. 

In general, at very early times, at very high temperatures, other particle 
degrees of freedom join the radiation background when T~mj for each 
particle type i, if that type is brought into thermal equilibrium through 
interactions. In equilibrium the energy density of a particle type i is given 
by 

Pi = j Eidn qi (12) 

and 

dn qi = ^[exp[(E qi - ^)/T] ± l]-\ 2 dq (13) 
where again gi counts the total number of degrees of freedom for type i, 

E qi = (m 2 i + q ff' 2 (14) 

Hi is the chemical potential if present and ± corresponds to either Fermi 
or Bose statistics. 

In the limit that T 3> rrii the total energy density can be conveniently 
expressed by 

where gs(F) are the total number of boson (fermion) degrees of freedom and 
the sum runs over all boson (fermion) states with m ^ T. The factor of 
7/8 is due to the difference between the Fermi and Bose integrals. Equation 
(15) defines N(T) by taking into account new particle degrees of freedom 
as the temperature is raised. 

In the radiation dominated epoch, we can obtain a relationship between 
the age of the Universe and its temperature 

( 90 y/ 2 2 

Put into a more convenient form 

tT 2 MeV = 2A[N(T)]- 1 / 2 (17) 
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TABLE 1. Effective numbers of degrees of freedom in the standard 
model. 



Temperature 


New Particles 


4N(T) 


T < m e 


•y's + i^'s 


29 


m e <T < m M 


e± 


43 


m M <T < 




57 


<T < Tc* 


7r's 


69 


T c <C T <C 7Tt s trangc 


- 7r's + U, M, d, ( 


i + gluons 205 


TYls <C T <C TYlcharm 


S, S 


247 


m c < T < m T 


c, c 


289 


m T < T < rribottom 


r± 


303 


nib <T < mw,z 


b,b 


345 


m w ,z <T < rritop 


W ± ,Z 


381 


m t <T < TUHiggs 


t,t 


423 




H° 


427 



*T C corresponds to the confmement-deconnnement transition between quarks and 
hadrons. N(T) is shown in Figure 1 for T c — 150 and 400 MeV. It has been assumed 
that niHiggs > rritop- 

where t is measured in seconds and TMeV m units of MeV. 

The value of N(T) at any given temperature depends on the particle 
physics model. In the standard 51/(3) x SU(2) x U (1) model, we can specify 
N(T) up to temperatures of 0(100) GeV. The change in N can be seen in 
the following table. 

At higher temperatures (T » 100 GeV), N(T) will be model dependent. 
For example, in the minimal SU(5) model, one needs to add to N(T), 24 
states for the X and Y gauge bosons, another 24 from the adjoint Higgs, 
and another 6 (in addition to the 4 already counted in W ± , Z and H) from 
the 5 of Higgs. Hence for T > M x in minimal SU(5), N(T) = 160.75. 
In a supersymmetric model this would at least double, with some changes 
possibly necessary in the table if the lightest supersymmetric particle has 
a mass below M#. 

The presence of a particle species in the thermal background assumes 
thermal equilibrium and hence interaction rates which are sufficiently fast 
compared with the expansion rate of the Universe. Roughly, this translates 
to the condition for each particle type i, that some rate Tj involving that 
type be larger than the expansion rate of the Universe or 

T l >H (18) 

in order to be in thermal equilibrium. 
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1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 

Log(TYMeV) 

Figure 1. The effective numbers of relativistic degrees of freedom as a function of 
temperature. 



Examples of a processes in equilibrium at early times which drops out 
of equilibrium or decouples at later times are the processes which involve 
neutrinos. If we consider the standard neutral or charged-current interac- 
tions such as e + + e~ ^v + voie + v^e + v etc., the rates for these 
processes can be approximated by 

r = n(av) (19) 

where (av) is the thermally averaged weak interaction cross section 

(<tv)~ 0(1(T 2 )T 2 /M^ (20) 

and n is the number density of leptons. Hence the rate for these interactions 
is 

T wk ~ 0(1(T 2 )T 5 /M^ (21) 
The expansion rate, on the other hand, is just 

\ 1/2 



H = (^Jlpy 2 = (^-N(T) j T 2 /M P ~ l.mN{TY/ 2 T 2 /M P . 



The Planck mass M P = G 



-1/2 



N 



1.22 x 10 19 GeV. 



Neutrinos will be in equilibrium when T w k > H or 
T > (bmM^/Mp) 1 ^^ IMeV. 



(22) 



(23) 
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The temperature at which these rates are equal is commonly referred to as 
the decoupling or freeze-out temperature and is defined by 

T(T d ) = H(T d ) (24) 

For temperatures T > T d , neutrinos will be in equilibrium, while for T <T d 
they will not. Basically, in terms of their interactions, the expansion rate 
is just too fast and they never "see" the rest of the matter in the Universe 
(or themselves). Their momenta will simply redshift and their effective 
temperature (the shape of their momenta distribution is not changed from 
that of a blackbody) will simply fall with T~l/R. 

Soon after decoupling the e^ pairs in the thermal background begin 
to annihilate (when T < m e ). Because the neutrinos are decoupled, the 
energy released heats up the photon background relative to the neutri- 
nos. The change in the photon temperature can be easily computed from 
entropy conservation. The neutrino entropy must be conserved separately 
from the entropy of interacting particles. If we denote T>, the temperature 
of photons, and before annihilation, we also have T u = T> as well. The 
entropy density of the interacting particles at T = T> is just 

while at T = T< , the temperature of the photons just after e ± annihilation, 
the entropy density is 

,26, 

and by conservation of entropy s < = s > and 

(T</T>) 3 = 11/4 (27) 

Thus, the photon background is at higher temperature than the neutrinos 
because the annihilation energy could not be shared among the neutri- 
nos, and 

T u = (4/H) 1 / 3 ^ ~ 1.9K (28) 
1.2. HISTORICAL PERSPECTIVES 

There has always been an intimate connection between BBN and the mi- 
crowave background as a key test to the standard big bang model. Indeed, 
it was the formulation of BBN which predicted the existence of the mi- 
crowave background radiation [4]. The argument is rather simple. BBN 
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requires temperatures greater than 100 keV, which according to eqs. (16) 
and (17) corresponds to timescales less than about 200 s. The typical cross 
section for the first link in the nucleosynthetic chain is 



av(p + D + 7) ~ 5 x 10 20 cm 3 /s 
This implies that it was necessary to achieve a density 



n 



— ~ 10 17 cm- 3 



avt 



(29) 



(30) 



The density in baryons today is known approximately from the density of 
visible matter to be ng ~ 10~ 7 cm" 3 and since we know that that the 
density n scales as R~ 3 ~ T 3 , the temperature today must be 



T = (n Bo /n) 1 / 3 T BBN ~ 10K 



(31) 



A pretty good estimate. 




Figure 2. The helium (Y) vs oxygen (O/H) abundances in extragalactic HII regions 
emphasizing the lack of low 4 He regions. 

Despite its simplicity, BBN was criticized early on, due to its shortcom- 
ings in being able to produce the observed abundances of all of the element 
isotopes due primarily to the gaps in stable nuclei at A = 5 and A = 8. At- 
tention was therefore turned to stellar nucleosynthesis [5]. However, while 
the elements from helium on up can be and are produced in stars, no other 
astrophysical site has ever survived for the production of deuterium [6] . In 
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addition, if one assumes that 4 He is entirely of stellar origin, one should 
be able to find places in the Universe in which the 4 He mass fraction is 
substantially below 25%. The 4 He data shown in Figure 2, emphasizes the 
fact that indeed no such region with low 4 He has ever been observed and 
that (together with the need to produce D) leads one to conclude that 
BBN nucleosynthesis is a necessary component in any cosmological model. 
The foundations of modern BBN continued to be laid over time [7] estab- 
lishing the notions of equilibrium and the nuclear network to obtain the 
abundances of D through 7 Li. 

2. Theory 

Conditions for the synthesis of the light elements were attained in the early 
Universe at temperatures T <> 1 MeV. At somewhat higher temperatures, 
weak interaction rates were in equilibrium. In particular, the processes 

p + v e 
p + e~ 
V + e~ + v e 

fix the ratio of number densities of neutrons to protons. At T > 1 MeV, 
(n/p) ~ 1. The energy density (dominated by radiation) in the standard 
model is 

' = li ! V>' < 32 > 

from the contributions of photons, electrons and positrons, and N v neutrino 
flavors. 

As the temperature fell and approached the point where the weak in- 
teraction rates were no longer fast enough to maintain equilibrium, the 
neutron to proton ratio was given approximately by the Boltzmann factor, 
(n/p) ~ e -Am/T^ w h ere Am is the neutron-proton mass difference. The 
final abundance of 4 He is very sensitive to the (n/p) ratio. As in the case 
of the neutrino interactions discussed above, freeze out occurs at about an 
MeV (slightly less than an MeV in this case). 

The nucleosynthesis chain begins with the formation of deuterium through 
the process, p + n — > D +7. However, because the large number of photons 
relative to nucleons, n~ l = n^/ns ~ 10 10 , deuterium production is delayed 
past the point where the temperature has fallen below the deuterium bind- 
ing energy, Eb = 2.2 MeV (the average photon energy in a blackbody is 
E^ ~ 2.7T). The point being that there are many photons in the exponen- 
tial tail of the photon energy distribution with energies E > Eb despite 
the fact that the temperature or E 1 are less than Eb- This can be seen 



n + e <-> 
n + v e <-> 
n <-> 
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by comparing the qualitative expressions for the deuterium production and 
destruction rates, 

T p ps ub<jv (33) 
T d ps nj ave~ EB/T 

When the quantity 77 _1 exp(— Eg/T) ~ 1 the rate for deuterium destruction 
(D +7 — > p + n) finally falls below the deuterium production rate and the 
nuclear chain begins at a temperature T ~ O.lMeV. 

In addition to the p (n, 7) D reaction, the other major reactions leading 
to the production of the light elements are: 

D (D, p) T D (n, 7) T 3 He (n, p) T 

D (D, n) 3 He D (p, 7) 3 He 

Followed by the reactions producing 4 He: 

D (D, 7) 4 He 3 He ( 3 He, 2p) 4 He 
D ( 3 He, p) 4 He T (p, 7) 4 He 
T (D, n) 4 He 3 He (n, 7) 4 He 

The gap at A = 5 is overcome and the production of 7 Li proceeds through: 

3 He ( 4 He, 7 ) 7 Be 

-> 7 Li + e+ + v e 

T ( 4 He,7) 7 Li 

The gap at A = 8 prevents the production of other isotopes in any signifi- 
cant quantity The nuclear chain in BBN calculations was extended [8] and 
is shown in Figure 3. 

The dominant product of big bang nucleosynthesis is 4 He resulting in 
an abundance of close to 25% by mass. This quantity is easily estimated 
by counting the number of neutrons present when nucleosynthesis begins. 
When the weak interaction rates responsible for n — p equilibrium freeze- 
out, at T Ps 0.8 MeV, the neutron to proton ratio is about 1/6. When free 
neutron decays prior to deuterium formation are taken into account, the 
ratio drops to (n/p) ps 1/7. Then simple counting yields a 4 He mass fraction 

In the standard model with N u = 3, there is basically one free parameter 
in BBN, namely the baryon to photon ratio, r]. As we have seen above, 
the value of r] controls the onset of nucleosynthesis through the deuterium 
bottleneck. For larger values of rj, the quantity ?7~ 1 exp(— Eb/T) is smaller, 
and hence the nuclear chain may begin at a higher temperature. Remember 
also that a key ingredient in determining the final mass fraction of 4 He, is 
(n/p) [see eq. (34)] and that the final value of (n/p) was determined by the 



11 




Figure 3. The nuclear network used in BBN calculations. 



time at which nucleosynthesis begins, thus controlling the time available for 
free decays after freeze out. If nucleosynthesis begins earlier, this leaves less 
time for neutrons to decay and the value of (n/p) and hence Y p is increased. 
But because the (n/p) ratio is only weakly dependent on v, the 4 He mass 
fraction is relatively flat as a function of 77. When we go beyond the standard 
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Figure 4- The light element abundances from big bang nucleosynthesis as a function of 
r/, including 6 Li. 

model, the 4 He abundance is very sensitive to changes in the expansion 
rate which can be related to the effective number of neutrino flavors as 
will be discussed below. Lesser amounts of the other light elements are 
produced: D and 3 He at the level of about 10 -5 by number, and 7 Li at the 
level of 10 -10 by number. These abundances (along with 6 Li) are shown in 
Figure 4 [8]. In Figure 5, the produced abundances of the intermediate mass 
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v 



Figure 5. The intermediate mass element abundances from big bang nucleosynthesis as 
a function of rj. 



isotopes 9 Be, °B, n B are also shown. These abundances are far below the 
observed values and it is believed that these isotopes are formed in cosmic 
ray nucleosynthesis. 

It is perhaps convenient at this time to note that the value of 77 is directly 
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related to the fraction of f2 in baryons. Indeed, one can write 

n B h 2 = 3.67 x 10 7 r/(T /2.728K) 3 (35) 

where To is the present temperature of the microwave background. 

Historically, it has been common to refer to two other parameter in 
BBN, the neutron mean life and the number of neutrino flavors. The neu- 
tron mean life is now very well determined and its remaining uncertainty 
can be treated simply as an uncertainty in the calculated abundance of 4 He. 
Although the number of neutrino flavors has also been fixed experimentally, 
BBN is sensitive to the number of light degrees of freedom whether or not 
they interact weakly. It is often convenient to refer to these degrees of free- 
dom as neutrino equivalents. By increasing N v in eq. (32), one increases 
the expansion rate H oc ^fp. As a result, the weak interactions freeze-out 
at a higher temperature (see eqs. (21) and (22)). This leads once again to 
a higher value for (n/p) and a higher 4 He mass fraction. This effect will be 
treated in more detail in the last lecture. 

For the comparison with the observations, I will use the resulting abun- 
dances of the light elements shown in Figure 6, which concentrate on the 
range in rjiQ between 1 and 10. The curves for the He mass fraction, Y, 
bracket the computed range based on the uncertainty of the neutron mean- 
life which has been taken as [9] r ra = 887±2 s. Uncertainties in the produced 
7 Li abundances have been adopted from the results in Hata et al. [10]. Un- 
certainties in D and 3 He production are small on the scale of this figure. 
The boxes correspond to the observed abundances and will be discussed 
below. 

3. Abundances 

3.1. 4 Hc 

4 He is produced in stars, therefore to be able to extract a primordial 
abundance of 4 He, it is advantageous to make abundance measurements in 
systems of very low metallicity. Low metallicity or low abundances of C, 
N, and O relative to the solar abundances of these elements would indicate 
that the level of stellar processing is minimized. The 4 He abundance in very 
low metallicity regions is best determined from observations of Hell — > Hel 
recombination lines in extragalactic HII (ionized hydrogen) regions. There is 
now a good collection of abundance information on the 4 He mass fraction, 
Y, O/H, and N/H in over 70 [11, 12, 13] such regions. In an extensive 
study based on the data in [11, 12], it was found [14] that the data is well 
represented by a linear correlation for Y vs. O/H and Y vs. N/H. It is then 
expected that the primordial abundance of 4 He can be determined from 
the intercept of that relation. A detailed analysis of the data including 
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0.26 




^10 

Figure 6. The light element abundances from big bang nucleosynthesis as a function of 
7710 = 10 1C V 



that in [13] found an intercept corresponding to a primordial abundance 
Y p = 0.234 ± 0.002 ± 0.005 [15]. The stability of this fit was verified by 
a statistical bootstrap analysis [16] showing that the fits were not overly 
sensitive to any particular HII region. 

To make use of the 4 He data, it is crucial to obtain high quality and very 
low metallicity data. In principle, any one HII region (with non-zero metal- 
licity) should provide an upper limit to Y p since some stellar processing has 
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Figure 7. Equal likelihood 1, 2, and 3 a contours in the w — Y p plane. The cross shows 
the position of the peak of the likelihood function. 



taken place augmenting the primordial value. Thus the determination of 
Y p by an extrapolation to zero metallicity could be avoided by the obser- 
vations of either low metallicity or low helium HII regions. For a very low 
metallicity HII region such an upper limit may even provide a reasonable 
estimate of Y p . 

Another way to avoid an extrapolation to zero metallicity (though such 
an extrapolation is in fact quite minimal given the low metallicity data 
available), one can perform a Bayesian analysis [17] in which one makes 
no other assumption other than the observed 4 He abundance is greater 
than or equal to the primordial abundance. That is, one can consider three 
quantities: Yt, the true 4 He abundance in an H II region about which 
the observed abundance Yo is distributed. Both of these differ from the 
primordial abundance Y p and the only prior assumed is that for each object, 
Yr > Y p . If we assume that the true abundance differs from the primordial 
abundance by no more than w, we can derive a total likelihood function 
by integrating out the unknown Yp. We can then plot the equal likelihood 
contours as a function of Y p and w . This is shown in Figure 7, where the 
32 points of lowest metallicity have been used to calculate the likelihood 
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Figure 8. The helium (Y) and oxygen (O/H) abundances in extragalactic HII regions, 
from refs. [11, 12] (open circles), and from ref. [18] (filled circles). Lines connect the same 
regions observed by different groups. 

function. The peak occurs at Y p = 0.238 and the most likely width is 
w = 0.009. The 95% CL upper limit to Y p in this case is 0.245. For further 
details on this approach see [17]. The data used in Figure 7 represents an 
update of that work and includes the data of ref. [18]. 

Although the above estimates on Y p are consistent with those based 
on a linear extrapolation of the data, it has been claimed that the new 
data in refs. [13] and [18] leads to a significantly higher value for Y p (in 
excess of 24%). The higher values of Y p , quoted by Izotov and Thuan are 
based on their 4 He abundances derived by their method of determining all 
of the parameters from a set of 5 helium recombination lines. This gives 
Y p = 0.2444 ± 0.0015 + (44 ± 19)0/H. However, as argued in [15] there are 
inherent uncertainties in this method which are not reflected in the error 
budget. For this reason and because we can more easily compare their data 
with previous data, we use their results which are based on S II densities. 
These results are entirely consistent with the data in refs. [11, 12] as can 
be seen in Figure 8 where the Y versus O/H data from refs. [11, 12] (open 
circles) is shown along with the newer data (filled circles, from ref. [18]). 
The fit to the open circles is shown by the dashed line with intercept 0.234, 
the fit to the filled circles is shown by the this solid line (barely visible) 
with intercept 0.239. Combining all the data one finds [19] a 4 He mass 
fraction based on 73 distinct HII regions 



Y v = 0.238 ± 0.002 ± 0.005 



(36) 
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The first uncertainty is purely statistical and the second uncertainty is an 
estimate of the systematic uncertainty in the primordial abundance de- 
termination [15]. The fit to all the data is shown by the thick line with 
intercept given by Eq. (36) above. The small errors quoted in [18] account 
for the total fit being skewed to the higher value of Y p . The solid box for 
4 He in Figure 6 represents the range (at 2u s t a t) from (36). The dashed box 
extends this by including the systematic uncertainty. A somewhat lower 
primordial abundance of Y p = 0.235 ± .003 is found by restricting to the 36 
most metal poor regions [19]. 

The primordial 4 He abundance can also be determined by examining the 
correlation between Y and N/H. Indeed in all but one of the H II regions, 
N/H data is also available. However, unless N/H is directly proportional to 
O/H, it is not clear that a linear Y vs. N/H fit should give the same result. 
Indeed, the proportionality of N/H to O/H (or in other terms the primary vs 
secondary nature of nitrogen) has been studied [14, 15, 19]. Unfortunately 
from a theoretical point of view this question lies in the realm of very 
uncertain yields for nitrogen in AGB stars. The data indicate that N is 
mostly primary. Though the secondary contribution may be responsible for 
yielding systematically higher intercepts for Y vs N/H relative to Y vs. 
O/H, however the difference is small < 0.003. 

Finally, it also of interest to test our understanding of the slope in the 
Y vs. O/H data. The data overall show a relatively steep slope AY/ AO ~ 
110 ± 25. Models of chemical evolution typically give a much smaller value 
of about 20 and even in models with outflow (material ejected from the 
galaxy) the slopes only go up to about 60. This question, like the N vs. O 
question is highly sensitive to very uncertain theoretical yields [19]. 

3.2. 7 Li 

The 7 Li abundance is also reasonably well known. In old, hot, population- 
II stars, 7 Li is found to have a very nearly uniform abundance [20]. For 
stars with a surface temperature T > 5500 K and a metallicity less than 
about l/20th solar (so that effects such as stellar convection may not be 
important), the abundances show little or no dispersion beyond that which 
is consistent with the errors of individual measurements. Indeed, as detailed 
in ref. [21], much of the work concerning 7 Li has to do with the presence 
or absence of dispersion and whether or not there is in fact some tiny slope 
to a [Li] = log 7 Li/H + 12 vs. T or [Li] vs. [Fe/H] relationship ([Fe/H] is 
the log of the Fe/H ratio relative to the solar value). 

There is 7 Li data from nearly 100 halo stars, from a variety of sources. 
When the Li data from stars with [Fe/H] < -1.3 is plotted as a function of 
surface temperature, one sees a plateau emerging for T > 5500 K as shown 
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Figure 9. The Li abundance in halo stars with [Fe/H] < -1.3, as a function of surface 
temperature. The dashed line shows the value of the weighted mean of the plateau data. 

in Figure 9 for the data taken from ref. [21]. As one can see from the 
figure, at high temperatures, where the convection zone does not go deep 
below the surface, the Li abundance is uniform. At lower temperatures, the 
surface abundance of Li is depleted as Li passes through the hotter interior 
of the star and is destroyed. The lack of dispersion in the plateau region is 
evidence that this abundance is indeed primordial (or at least very close to 
it). Another way to see the plateau is to plot the Li abundance data as a 
function of metallicity, this time with the restriction that T > 5500 K as 
seen in Figure 10. In this case, the plateau emerges at low [Fe/H] as would 
be expected. At higher [Fe/H], the convection zone remains near the surface 
only for much hotter stars. Thus, for [Fe/H] > -1.3, the effects of depletion 
are seen. Also apparent in this figure is that at higher metallicities there is 
evidence for the production of Li which rises by over an order of magnitude 
at solar metallicity. 

I will use the value given in ref. [21] as the best estimate for the mean 
7 Li abundance and its statistical uncertainty in halo stars 

Li/H = (1.6 ± 0.1) x 10~ 10 (37) 

The Li abundance determination is sensitive to stellar parameters such as 
the assumed surface temperature, the metallicity and the surface gravity. 
The greatest model dependence is on the conversion of the observed colors 
(B-V) to temperature. For example, in a sample of 55 stars taken from 
the papers of ref. [20], one finds [Li] = 2.08 ± 0.02. From Thorburn's[22] 
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Figure 10. The Li abundance in halo stars with T > 5500 K, as a function of metallicity. 
Filled diamonds represent upper limits. 



sample of 74 stars one finds [Li] = 2.26 db 0.01. I believe that much of 
this systematic uncertainty has now been relieved by newer methods for 
determining the stellar parameters and the Li abundance. The two papers 
in ref. [21], are based on different methods. The first uses temperatures 
which are determined by observations of Balmer lines [23] and the second 
employs the infra-red flux method [24]. The data sets in these papers which 
contain 24 and 41 stars respectively, both give [Li] = 2.21 ± 0.01. The Li 
abundance used in (37) corresponds to this value. In fact when correlated 
with respect to either temperature or iron, the data in the latter paper of 
[21] shows no statistical trend. With respect to temperature, the data give 
[Li] = 1.28 ± 0.43 + (0.015 + 0.007)T/100K, indicating a change in [Li] well 
within the data over the observed temperature range. With respect to iron, 
no slope is found [Li] = 2.17 ± 0.07 + (-0.018 ± 0.031) [Fe/H]. The variance 
in the plateau data is less than 0.01. The solid box for 7 Li in Figure 6 
represents the 2<7 s t a t range from (37). 

There is however an important source of systematic error due to the 
possibility that Li has been depleted in these stars from their initial abun- 
dance. These uncertainties are however limited. As was noted above, the 
lack of dispersion in the Li data limits the amount of depletion. In addi- 
tion, standard stellar models [25] predict that any depletion of 7 Li would 
be accompanied by a very severe depletion of 6 Li. Until recently, 6 Li had 
never been observed in hot pop II stars. The observation [26] of 6 Li (which 
turns out to be consistent with its origin in cosmic-ray nucleosynthesis and 
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with a small amount of depletion as expected from standard stellar mod- 
els) is another good indication that 7 Li has not been destroyed in these 
stars[27, 28, 29]. 

Aside from the big bang, Li is produced together with Be and B in 
cosmic ray spallation of C,N,0 by protons and a-particles. Li is also pro- 
duced by a — a fusion. Be and B have been observed in these same pop II 
stars and in particular there are a dozen or so stars in which both Be and 
Li have been observed. Thus Be (and B though there is still a paucity of 
data) can be used as a consistency check on primordial Li [30]. Based on 
the Be abundance found in these stars, one can conclude that no more than 
10-20% of the 7 Li is due to cosmic ray nucleosynthesis leaving the remain- 
der (an abundance near 10 -10 ) as primordial. The dashed box in Figure 6, 
accounts for the possibility that as much as half of the primordial 7 Li has 
been destroyed in stars, and that as much as 20% of the observed 7 Li may 
have been produced in cosmic ray collisions rather than in the Big Bang. 
For 7 Li, the uncertainties are clearly dominated by systematic effects. 



3.3. D 

Turning to D/H, we have three basic types of abundance information: 1) 
ISM data, 2) solar system information, and perhaps 3) a primordial abun- 
dance from quasar absorption systems. The best measurement for ISM D/H 
is [31] 

(D/H)i SM = 1-60 ± 0.09lHo x !0~ 5 (38) 

Because there are no known astrophysical sites for the production of deu- 
terium, all observed D must be primordial. As a result, a firm lower bound 
from deuterium establishes an upper bound on rj which is robust and is 
shown by the lower right of the solid box in Figure 6. The solar abundance 
of D/H is inferred from two distinct measurements of 3 He. The solar wind 
measurements of 3 He as well as the low temperature components of step- 
wise heating measurements of 3 He in meteorites yield the presolar (D + 
3 He)/H ratio, as D was efficiently burned to 3 He in the Sun's pre-main- 
sequence phase. These measurements indicate that [32, 33] 

(rj i 3tt„\ 
j = (4.1 ±0.6 ±1.4) x 10" 5 (39) 

The high temperature components in meteorites are believed to yield the 
true solar 3 He/H ratio of [32, 33] 

/ 3 He\ 

"' ' = (1.5 ±0.2 ±0.3) x 10" 5 (40) 



H 

/ 
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The difference between these two abundances reveals the presolar D/H 
ratio, giving, 



It should be noted that measurements of surface abundances of HD on 
Jupiter show a somewhat higher value for D/H, D/H = 5 ± 2 x 10~ 5 [34]. 
If this value is confirmed and if fractionation does not significantly alter 
the D/H ratio (as it was suspected to for previous measurements involving 
CH3D), it may have an important impact on galactic chemical evolution 
models. This value is marginally consistent with the inferred meteoritic 
values. 

Finally, there have been several reported measurements of D /H in high 
redshift quasar absorption systems. Such measurements are in principle ca- 
pable of determining the primordial value for D/H and hence rj, because of 
the strong and monotonic dependence of D/H on rj. However, at present, 
detections of D/H using quasar absorption systems do not yield a conclu- 
sive value for D/H. As such, it should be cautioned that these values may 
not turn out to represent the true primordial value and it is very unlikely 
that both are primordial and indicate an inhomogeneity [35] (a large scale 
inhomogeneity of the magnitude required to placate all observations is ex- 
cluded by the isotropy of the microwave background radiation). The first 
of these measurements [36] indicated a rather high D/H ratio, D/H « 1.9 
- 2.5 xlO~ 4 . Other high D/H ratios were reported in [37]. More recently, a 
similarly high value of D/H = 2.0 ±0.5 x 10 -4 was reported in a relatively 
low redshift system (making it less suspect to interloper problems) [38]. 
However, there are reported low values of D/H in other such systems [39] 
with values of D/H originally reported as low as ~ 2.5 x 10~ 5 , significantly 
lower than the ones quoted above. The abundance in these systems has 
been revised upwards to about 3.4 ±0.3 x 10~ 5 [40]. I will not enter into 
the debate as to which if any of these observations may be a better repre- 
sentation of the true primordial D/H ratio. I only note that it remains a 
highly contested issue [40, 41] The range of quasar absorber D/H is shown 
by the dashed box in Figure 6. 

There are also several types of 3 He measurements. As noted above, 
meteoritic extractions yield a presolar value for 3 He/H as given in Eq. 
(40). In addition, there are several ISM measurements of 3 He in galactic 
HII regions [42] which show a wide dispersion which may be indicative of 
pollution or a bias [43] 



(D/H) « (2.6 ± 0.6 ± 1.4) x 10~ 5 



(41) 




(42) 
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There is also a recent ISM measurement of 3 He [44] with 

(^) = 2.H:I x 10- 5 (43) 
V H /ism 

Finally there are observations of 3 He in planetary nebulae [45] which show 
a very high 3 He abundance of 3 He/H ~ 1CT 3 . 

Each of the light element isotopes can be made consistent with theory 
for a specific range in r\. Overall consistency of course requires that the range 
in rj agree among all four light elements. However, as will be argued below 
D and 3 He are far more sensitive to chemical evolution than 4 He or 7 Li and 
as such the direct comparison between the theoretical predictions of the 
primordial abundances of D and 3 He with the observational determination 
of their abundances is far more difficult. Therefore in what follows I will for 
the most part restrict the comparison between theory and observation to 
the two isotopes who suffer the least from the effects of chemical evolution. 



4. Chemical Evolution 

Because we can not directly measure the primordial abundances of any of 
the light element isotopes, we are required to make some assumptions con- 
cerning the evolution of these isotopes. As has been discussed above, 4 He 
is produced in stars along with oxygen and nitrogen. 7 Li can be destroyed 
in stars and produced in several (though still uncertain) environments. D 
is totally destroyed in the star formation process and 3 He is both produced 
and destroyed in stars with fairly uncertain yields. It is therefore preferable, 
if possible to observe the light element isotopes in a low metallicity envi- 
ronment. Such is the case with He and 7 Li, and we can be fairly assured 
that the abundance determinations of these isotopes are close to primordial. 
If the quasar absorption system measurements of D/H stabilize, then this 
too may be very close to a primordial measurement. Otherwise, to match 
the solar and present abundances of D and 3 He to their primordial values 
requires a model of galactic chemical evolution. 

The main inputs to chemical evolution models are: 1) The initial mass 
function, <f)(m), indicating the distribution of stellar masses. Typically, a 
simple power law form for the IMF is chosen, 4>{m) ~ m~ x , with x ~ —2.7. 
This is a fairly good representation of the observed distribution, particularly 
at larger masses. 2) The star formation rate, ip. Typical choices for a SFR 
are ip(t) oc a or a 2 or even a straight exponential e~*/ r . a is the fraction 
of mass in gas, M gas /M tot . 3) The presence of infalling or outflowing gas; 
and of course 4) the stellar yields. It is the latter, particularly in the case of 
3 He, that is the cause for so much uncertainty. Chemical evolution models 
simply set up a series of evolution equations which trace desired quantities. 



24 



For example, the mass in gas and the SFR evolve through a relation such 
as 

^ = -m + e(t)+i(t)-o(t) (44) 

where e represents the amount of gas ejected from stars, i is the gas infall 
rate, and o is the gas outflow rate. The ejection rate is in turn given by 

e(t) = J (m — mji)4>(m)ip(t — T(m))dm (45) 

where ran is the remnant mass (a function of the stellar mass m as well) 
and r(m) is the stellar lifetime. If we ignore r(m), then the ejection rate is 
simply proportional to the star formation rate ip, e{t) = Rip- R is referred to 
as the return fraction and this approximation is known as the instantaneous 
recycling approximation (IRA). Similar equations can be developed which 
trace the abundances of the element isotopes [46]. Neglecting both infall 
and outflow, these take the form 



d{XM gas ) 
dt 



-^{t)X + e x {t) (46) 



where X is the mass fraction of a particular element of interest and ex is 
the mass fraction of the element ejected in the death of a star. In the case 
of deuterium, = 0. 

As one can see from (46) deuterium is always a monotonically decreasing 
function of time in chemical evolution models. The degree to which D is 
destroyed, is however a model dependent question which depends sensitively 
on the IMF and SFR. The evolution of 3 He is however considerably more 
complicated. Stellar models predict that substantial amounts of 3 He are 
produced in stars between 1 and 3 M . For M < 8M , Iben and Truran 
[47] calculate 

( 3 He/H) / = 1.8 x 10~ 4 {^^j 2 + 0.7 [(D + 3 He) /h] . (47) 

so that for example, when r/io = 3, ((D + 3 He)/H)j = 9 x 10~ 5 , and the 
ratio of the final abundance of 3 He/H to the initial (D + 3 He) /H abundance 
denoted by 53 is gs(lM Q ) = 2.7. The 3 He abundance is nearly tripled. It 
should be emphasized that this prediction is in fact consistent with the 
observation of high 3 He/H in planetary nebulae [45]. 

Generally, implementation of the 3 He yield in Eq. (47) in chemical evo- 
lution models leads to an overproduction of 3 He/H particularly at the solar 
epoch [43, 48]. For example, in Figure 11, the evolution of D and 3 He is 
shown for a model in which only a modest amount of deuterium is de- 
stroyed. Namely, by a factor of 5, from D/H = 7.5 xl0~ 5 to a present 
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Figure 11. The evolution of D/H and 3 He/H and (D+ 3 He)/H with time in units of 
10~ 5 . The assumed primordial abundance of D/H is 7.5 xlO -5 . The solid curve shows 
the evolution of 3 He/H, the dashed curve for D/H and the dotted curve for the sum 
(D+ 3 He)/H. The diamonds show the data for 3 He, the open squares for deuterium and 
the open circle for the sum. 



value of ~ 1.5 x 10 -5 . However, due to the production of 3 He in low mass 
stars, 3 He is greatly overproduced relative to the solar value. This problem 
is compounded in models with an intense period of D destruction. In Scully 
et al. [49], a dynamically generated supernovae wind model was coupled to 
models of galactic chemical evolution with the aim of reducing a primordial 
D/H abundance of 2 xl0~ 4 to the present ISM value without overproduc- 
ing heavy elements and remaining consistent with the other observational 
constraints typically imposed on such models. In Figure 12, the evolution 
of D/H and 3 He/H is shown as a function of time in several representa- 
tive models with significant deuterium destruction factors (see ref [49] for 
details). However, as one can plainly see, 3 He is grossly overproduced (the 
deuterium data is represented by squares and 3 He by circles). 

The overproduction of 3 He relative to the solar meteoritic value seems 
to be a generic feature of chemical evolution models when 3 He production 
in low mass stars is included. This result appears to be independent of 
the chemical evolution model and is directly related to the assumed stellar 
yields of 3 He. It has recently been suggested that at least some low mass 
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stars may indeed be net destroyers of 3 He if one includes the effects of 
extra mixing below the conventional convection zone in low mass stars on 
the red giant branch [50, 51]. The extra mixing does not take place for stars 
which do not undergo a helium core flash (i.e. stars > 1.7 - 2 M ). Thus 
stars with masses less than 1.7 M Q are responsible for the 3 He destruction. 
Using the yields of Boothroyd and Malaney [51], it was shown [52] that 
these reduced 3 He yields in low mass stars can account for the relatively 
low solar and present day 3 He/H abundances observed. In fact, in some 
cases, 3 He was underproduced. To account for the 3 He evolution and the 
fact that some low mass stars must be producers of 3 He as indicated by the 
planetary nebulae data, it was suggested that the new yields apply only 
to a fraction (albeit large) of low mass stars [52, 53]. The corresponding 
evolution [52] of D/H and 3 He/H is shown in Figure 13. 

The models of chemical evolution discussed above indicate that it is 
possible to destroy significant amounts of deuterium and remain consistent 
with chemical evolutionary constraints. To do so however, comes with a 
price. Large deuterium destruction factors require substantial amounts of 
stellar processing, which at the same time produce heavy elements. To 
keep the heavy element abundances in the Galaxy in check, significant 
Galactic winds enriched in heavy elements must be incorporated. In fact 
there is some evidence that enriched winds were operative in the early 
Galaxy. In the X-ray cluster satellites observed by Mushotzky et al. [54] and 
Loewenstein and Mushotzky [55] the mean oxygen abundance was found to 
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be roughly half solar. This corresponds to a near solar abundance of heavy 
elements in the inter-Galactic medium, where apparently little or no star 
formation has taken place. 

If our Galaxy is typical in the Universe, then the models of the type dis- 
cussed above would indicate that the luminosity density of the Universe at 
high redshift should also be substantial augmented relative to the present. 
Recent observations of the luminosity density at high redshift [56] are mak- 
ing it possible for the first time to test models of cosmic chemical evolution. 
The high redshift observations, are very discriminatory with respect to a 
given SFR [57]. Models in which the star formation rate is proportional to 
the gas mass fraction (these are common place in Galactic chemical evo- 
lution) have difficulties to fit the multi-color data from z = to 1. This 
includes many of the successful Galactic infall models. In contrast, mod- 
els with a steeply decreasing SFR are favored. In Figure 14, the predicted 
luminosity density based on the model with evolution shown in Figure 13 
from [49], as compared with the observations (see ref. [57] for details). 

While it would be premature to conclude that all models with large 
deuterium destruction factors are favored, it does seem that models which 
do fit the high redshift data destroy significant amounts of D/H. On the 
other hand, we can not exclude models which destroy only a small amount 
of D/H as Galactic models of chemical evolution. In this case, however 
the evolution of our Galaxy is anomalous with respect to the cosmic av- 
erage. If the low D/H measurements [39, 40] hold up, then it would seem 



28 




Figure 14- The tricolor luminosity densities (UV, B and IR) at A = 0.28,0.44 and 1.0 
fim, m units of (h/.5) WHz^Mpc" 3 as a function of redshift for a model shown in 13 
which destroys significant amounts of D/H. The data are taken from [56]. 



that our Galaxy also has an anomalously high D/H abundance. That is 
we would predict in this case that the present cosmic abundance of D/H 
is significantly lower than the observed ISM value. If the high D/H obser- 
vations [36, 37, 38] hold up, we could conclude that our Galaxy is indeed 
representative of the cosmic star formation history. 
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5. Likelihood Analyses 

Monte Carlo techniques have proven to be a useful form of analysis for big 
bang nucleosynthesis [58, 59, 10]. An analysis of this sort was performed [2] 
using only 4 He and 7 Li. Two elements are sufficient for not only constraining 
the one parameter (77) theory of BBN, but also for testing for consistency. 
The procedure begins by establishing likelihood functions for the theory 
and observations. For example, for 4 He, the theoretical likelihood function 
takes the form 

^bbn(^) Ibbn) = e-^- y BB N (,)) 2 /2 CT ? (48) 

where Ibbn(^) is the central value for the 4 He mass fraction produced in 
the big bang as predicted by the theory at a given value of 77. u\ is the 
uncertainty in that value derived from the Monte Carlo calculations [10] 
and is a measure of the theoretical uncertainty in the big bang calculation. 
Similarly one can write down an expression for the observational likelihood 
function. Assuming Gaussian errors, the likelihood function for the obser- 
vations would take a form similar to that in (48). 

A total likelihood function for each value of 77 is derived by convolving 
the theoretical and observational distributions, which for 4 He is given by 

i 4He totai(r?) = J dYL BBN (Y, Ybbn (v)) Lo(Y, Y q ) (49) 

An analogous calculation is performed [2] for 7 Li. The resulting likelihood 
functions from the observed abundances given in Eqs. (36) and (37) is 
shown in Figure 15. As one can see there is very good agreement between 
4 He and 7 Li in the range of 7710 ~ 1.5 - 5.0. The double peaked nature 
of the 7 Li likelihood function is due to the presence of a minimum in the 
predicted lithium abundance. For a given observed value of 7 Li, there are 
two likely values of 77. 

The combined likelihood, for fitting both elements simultaneously, is 
given by the product of the two functions in Figure 15 and is shown in 
Figure 16. The 95% CL region covers the range 1.55 < 7710 < 4.45, with the 
two peaks occurring at 7710 = 1.9 and 3.5. This range corresponds to values 

between 

0.006 < n B h 2 < .016 (50) 

For the lower value of Y p = 0.234 ± 0.002 ± 0.005 as quoted in [15], 
the 4 He peak is shifted to slightly lower values of 77 and sits on top of the 
low- 77 peak as shown in Figure 17. (The difference in the 7 Li likelihood 
distribution is due to the assumed uncertainty in the 7 Li abundance which 
is slightly higher than that in Figure 15.) The combined likelihood in this 
case is shown in Figure 18. From Figure 18 it is clear that 4 He overlaps the 
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Figure 15. Likelihood distribution for each of 4 He and 7 Li, shown as a function of r\. The 
one-peak structure of the 4 He curve corresponds to the monotonic increase of Y v with 77, 
while the two peaks for 7 Li arise from the minimum in the 7 Li abundance prediction. 
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Figure 16. Combined likelihood for simultaneously fitting 4 He and 7 Li, as a function of 



31 




lower (in 77) 7 Li peak, and so one expects that there will be concordance in 
an allowed range of rj given by the overlap region. This is what one finds 
in Figure 18, which does show concordance and gives a preferred value for 
ri, n 10 = 1.8± 2 4 4 corresponding to (at 95% CL) tt B h 2 = .007±;gg|. 

Thus, we can conclude that the abundances of 4 He and 7 Li are con- 
sistent, and select an 7710 range which overlaps with (at the 95% CL) the 
longstanding favorite range around 7710 = 3. Furthermore, by finding con- 
cordance using only 4 He and 7 Li, we deduce that if there is problem with 
BBN, it must arise from D and 3 He and is thus tied to chemical evolution 
or the stellar evolution of 3 He. The most model-independent conclusion is 
that standard BBN with N v = 3 is not in jeopardy. 

It is interesting to compare the results from the likelihood function of 
4 He and 7 Li with that of D/H. Since D and 3 He are monotonic functions of 
77, a prediction for 77, based on 4 He and 7 Li, can be turned into a prediction 
for D and 3 He. The corresponding 95% CL ranges are D/H = (4.3 — 25) x 
1(T 5 and 3 He/H = (1.2 - 2.6) x 1(T 5 . If we did have full confidence in 
the measured value of D/H in quasar absorption systems, then we could 
perform the same statistical analysis using 4 He, 7 Li, and D. To include 
D/H, one would proceed in much the same way as with the other two light 
elements. We compute likelihood functions for the BBN predictions as in 
Eq. (48) and the likelihood function for the observations. These are then 
convolved as in Eq. (49). Using D/H = (2.0 ± 0.5) x 10~ 4 as indicated in 
the high D/H systems, we can plot the three likelihood functions including 
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Figure 18. Combined likelihood for simultaneously fitting 4 He and 7 Li, as a function of 
77 from from Figure 17. 

L^ tal (r?) in Figure 19. It is indeed startling how the three peaks, for D, 
4 He and 7 Li are in excellent agreement with each other. In Figure 20, the 
combined distribution is shown. We now have a very clean distribution and 
prediction for 77, 7710 = l-Sl 1 ^ 6 corresponding to Qsh 2 = .007l;ogf. The 
absence of any overlap with the high-77 peak of the 7 Li distribution has 
considerably lowered the upper limit to rj. Overall, the concordance limits 
in this case are dominated by the deuterium likelihood function. 

If instead, we assume that the low value [40] of D/H = (3.4 ±0.3) x 10~ 5 
is the primordial abundance, then we can again compare the likelihood 
distributions as in Figure 19, now substituting the low D/H value. As one 
can see from Figure 21, there is now hardly any overlap between the D 
and the 7 Li and 4 He distributions. The combined distribution shown in 
Figure 22 is compared with that in Figure 20. Though one can not use this 
likelihood analysis to prove the correctness of the high D /H measurements 
or the incorrectness of the low D/H measurements, the analysis clearly 
shows the difference in compatibility between the two values of D/H and 
the observational determinations of 4 He and 7 Li. To make the low D/H 
measurement compatible, one would have to argue for a shift upwards in 
4 He to a primordial value of 0.247 (a shift by 0.009) which is not warranted 
at this time by the data, and a 7 Li depletion factor of about 2, which is 
close to recent upper limits to the amount of depletion [60, 28]. 

The implications of the resulting predictions from big bang nucleosyn- 
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Figure 19. As in Figure 15, with the addition of the likelihood distribution for D/H 
assuming "high" D/H. 
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Figure 20. Combined likelihood for simultaneously fitting 4 He and 7 Li, and D as a 
function of r/ from Figure 19. 
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Figure 21. As in Figure 19, with the likelihood distribution for low D/H. 
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Figure 22. Combined likelihood for simultaneously fitting 4 He and 7 Li, and low D/H as 
a function of r\. The dashed curve represents the combined distribution shown in Figure 
20. 



thesis on dark matter are clear. First, if = 1 (as predicted by inflation), 
and £Ib ^0.1 which is certainly a robust conclusion based on D/H, then 
non-baryonic dark matter is a necessity. Second, on the scale of small groups 



35 



of galaxies, which are expected to sample the dark matter in galactic halos, 
> 0.05. This value can be compared with the best estimate for Qb from 
equation (50). The low r\ peak in Figure 16 corresponds to Ugh 2 = 0.007 
which for h = 1/2 gives Qb = 0.028. In this event, some non-baryonic dark 
matter in galactic halos is required. On the other hand, the high r\ peak 
in 16 corresponds to O^/i 2 = 0.013 which for h = 1/2 gives £Ib = 0.051 
and may be consistent with halo densities (but only just). If we include 
the data from the high D/H measurements in QSO absorbers this conclu- 
sion is unchanged, only the low rj peak survives, and some non-baryonic 
dark matter is needed in galactic halos. In contrast [61], the low D/H mea- 
surements would imply that Og/j 2 = 0.019 allowing for the possibility that 
SIb ~ 0.08. In this case, no non-baryonic dark matter is required in galactic 
halos. However, I remind the reader that the low D/H value of 3.4 xl0~ 5 
is at present only barely consistent with either the observations of 4 He or 
7 Li and their interpretations as being primordial abundances. 

6. Constraints from BBN 

Limits on particle physics beyond the standard model are mostly sensitive 
to the bounds imposed on the 4 He abundance. As is well known, the 4 He 
abundance is predominantly determined by the neutron-to-proton ratio just 
prior to nucleosynthesis and is easily estimated assuming that all neutrons 
are incorporated into 4 He (see Eq. (34)). As discussed earlier, the neutron- 
to-proton ratio is fixed by its equilibrium value at the freeze-out of the weak 
interaction rates at a temperature Tf ~ 1 MeV modulo the occasional free 
neutron decay. Furthermore, freeze-out is determined by the competition 
between the weak interaction rates and the expansion rate of the Universe 

G F 2 T f 5 ~ r wk (T» = H(T f ) ~ ^G^NTf 2 (51) 

where ./V counts the total (equivalent) number of relativistic particle species. 
The presence of additional neutrino flavors (or any other relativistic species) 
at the time of nucleosynthesis increases the overall energy density of the 
Universe and hence the expansion rate leading to a larger value of Tf, (n/p), 
and ultimately Y p . Because of the form of Eq. (51) it is clear that just as 
one can place limits [62] on N, any changes in the weak or gravitational 
coupling constants can be similarly constrained (for a recent discussion see 
ref. [63]). In concluding this lecture, I will discuss the current constraint on 
N v the number of particle species (in neutrino units) and the limit on the 
strength of new interactions, if 3 right-handed (nearly massless) neutrinos 
are assumed to exist. 

In the past, 3 He (together with D) has stood out in its importance for 
BBN, because it provided a (relatively large) lower limit for the baryon- 
to-photon ratio [64], 7710 > 2.8. This limit for a long time was seen to be 
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essential because it provided the only means for bounding rj from below and 
in effect allows one to set an upper limit on the number of neutrino flavors 
[62], N u , as well as other constraints on particle physics properties. That is, 
the upper bound to N v is strongly dependent on the lower bound to rj. This 



is easy to see: given an observed value for Y p , for lower rj, the 4 He abundance 
drops, allowing for a larger N v , which would raise the 4 He abundance back 
to the observed value. However, for rj < 4 x 10 -11 , corresponding to Sl^/i 2 <> 
.001 — .002, which is not too different from galactic mass densities, there 
is no bound whatsoever on N v [65]. Of course, with the improved data on 
7 Li, we do have lower bounds on r\ which exceed 10 -10 . In fact, despite the 
uncertainty in the D/H abundance in quasar absorption systems, the high 
D/H values can certainly be regarded as an upper limit to primordial D/H, 
which also yield a lower limit to rf. 

Because, of new observations of D and 3 He, and the new theoretical 
work on chemical evolution sparked by these observations, the bound on 
N v which is tied directly to these isotopes, should be called into question. As 
described earlier, the limits due to 3 He are ultimately tied to the assumed 
yields of low mass stars. Using the reduced yields as depicted in Figure 
13, consistent values of 77 < 2.8 are certainly possible. Ultimately, as I have 
said repeatedly, D/H measurements in quasar absorption systems may soon 
resolve this issue. However, the lower values of 77, relax the bounds on the 
number of neutrino flavors. 

As discussed above, the limit on N u comes about via the change in the 
expansion rate given by the Hubble parameter, 

ff2= 8^ = ^ sM + T AJVJrl (52) 

when compared to the weak interaction rates. Here Asm refers to the stan- 
dard model value for N. At T ~ 1 MeV, A S m = 43/4. Additional degrees 
of freedom will lead to an increase in the freeze-out temperature eventu- 
ally leading to a higher 4 He abundance. In fact, one can parameterize the 
dependence of Y on N u by 

Y = 0.2262 + 0.0131(JV„ - 3) + 0.0135 In r] W (53) 

in the vicinity of 7710 ~ 2. Eq. (53) also shows the weak (log) dependence 
on 77. However, rather than use (53) to obtain a limit, it is preferable to use 
the likelihood method. 

Just as 4 He and 7 Li were sufficient to determine a value for 77, a limit on 
N v can be obtained as well [2, 66]. The likelihood approach utilized above 
can be extended to include N u as a free parameter. Since the light element 
abundances can be computed as functions of both rj and N u , the likelihood 
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function can be denned by [66] 

L BB MY, Ibbn) = e -0-- W,,iM) W (54) 

and 

^ Hc totai(r/, K) = J (IYLbbn (Y, y B BN (v, N„)) L Q (Y, Y ) (55) 

Again, similar expressions are needed for 7 Li and D. A three-dimensional 
view of the combined likelihood functions [66] is shown in Figure 23 which is 
based on the slightly lower value of Y p as used in Figure 17. For updated (but 
similar) results see [67] . In this case the high and low 77 maxima of Figure 
18, show up as peaks in the L — 77 — N v space (L47 when D/H is neglected 
and L247 when high D/H is included). The peaks of the distribution as 
well as the allowed ranges of 17 and N v are more easily discerned in the 
contour plot of Figure 24 which shows the 50%, 68% and 95% confidence 
level contours in the two likelihood functions. The crosses show the location 
of the peaks of the likelihood functions. L47 peaks at N u = 3.0, 7710 = 1.8 
(in agreement with our previous results [2]) and at N v = 2.3, 7710 = 3.6. 
The 95% confidence level allows the following ranges in 77 and N v 

1.6 < N v < 4.0 1.3 < 7710 < 5.0 (56) 

Note however that the ranges in 77 and N u are strongly correlated as is 
evident in Figure 24. Since the deuterium likelihood function picks out a 
small range of values of 77, largely independent of N u , its effect on L247 
is to eliminate one of the two peaks in L47. L247 also peaks at N v = 3.0, 
7710 = 1.8. In this case the 95% contour gives the ranges 

2.0 <N V < 4.1 1.4 < 7710 < 2.6 (57) 

Finally, in Figure 25, the resulting 95 % CL upper limit to N v is shown 
as a function of Yp for several different choices for the primordial value of 
7 Li/H [67]. 

One should recall that the limit derived above is not meant for neutrinos 
in the strictest sense. That is, the limit is only useful when applied to 
additional particle degrees of freedom which necessarily do not couple to 
the Z°. For very weakly interacting particles, one must take into account 
the reduced abundance of these particles at the time of nucleosynthesis [68]. 
As discussed in the first lecture, the number of neutrinos today is reduced 
relative to the number of photons by (T^/T 7 ) 3 = 4/11. For some new 
particle, %, which decoupled at > 1 MeV, the same argument based on 
the conservation of entropy tells us that 



(58) 
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Figure 23. The combined two-dimensional likelihood functions for simultaneously fitting 
4 He and 7 Li in the top panel, and including D in the lower one as functions of both rj 
and N v . 



Thus we can translate the bound on N u , which is really a bound on the 
additional energy density at nucleosynthesis 
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Figure 24- The top panel shows contours in the combined likelihood function for 4 He 
and 7 Li. The contours represent 50% (innermost), 68% and 95% (outermost) confidence 
levels. The crosses mark the points of maximum likelihood. Also shown is the equivalent 
result when D is included. 



for additional boson states with gs degrees of freedom and fermion states 
with gp degrees of freedom. At nucleosynthesis T = T v = T~ and the limit 
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3 L_, — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — ,_l 
0.225 0.230 0.235 0.240 0.245 0.250 



Figure 25. The 95 % CL upper limit to N v as a function of Y p for (Li/H) p = 1.6, 2.0, 
2.6, 3.2, and 4.1 xlO" 10 . 



N v < 4 becomes 

?Ef(^) 4 + Ef(^) 4 <i (60) 

Such a limit would allow a single additional scalar degree of freedom (which 
counts as |) such as a majoron. On the other hand, in models with right- 
handed interactions, and three right-handed neutrinos, the constraint is se- 
vere. The right-handed states must have decoupled early enough to ensure 
(T Ur /T Ul ) 4 < 1/3. The temperature of a decoupled state is easily deter- 
mined from (58). Three right-handed neutrinos would require N(Td) > 25, 
which from Figure 1 implies that > 140 MeV, conservatively assuming a 
QCD transition temperature of 150 MeV. If right-handed interactions are 
mediated by additional gauge interactions, associated with some scale Mz>, 
and if we assume that the right handed interactions scale as M z , with a 
standard model-like coupling, then the decoupling temperature of the right 
handed interactions is related to Mz> by 

i^f ~ (^) 4 (61) 

J-dL M Z 

which for T^l ~ 3 MeV ( a more accurate value that the 1 MeV estimate) 
and ~ 140 MeV, we find that the associated mass scale becomes M Z / > 
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Figure 26. Limits on neutrino-like degrees of freedom. 



1.6 TeV! In general this constraint is very sensitive to the BBN limit on 
N v . In Figure 26, the allowed number of neutrino degrees of freedom are 
shown as a function of their decoupling temperature for the case of N v < 4 
and N v < 3.3, shown for comparison. 

7. Summary 

To summarize on the subject of big bang nucleosynthesis, I would assert 
that one can conclude that the present data on the abundances of the 
light element isotopes are consistent with the standard model of big bang 
nucleosynthesis. Using the the isotopes with the best data, 4 He and 7 Li, it 
is possible to constrain the theory and obtain a best set of values for the 
baryon-to-photon ratio of 7710 and the corresponding value for 

1.55 < 7710 < 4.45 95%CL , , 

.006 < n B h 2 < .016 95%CL { ' 

For 0.4 < h < 1, we have a range .006 < £Ib < -10. This is a rather low 
value for the baryon density and would suggest that much of the galactic 
dark matter is non-baryonic [69]. These predictions are in addition consis- 
tent with recent observations of D/H in quasar absorption systems which 
show a high value. Difficulty remains however, in matching the solar 3 He 
abundance, suggesting a problem with our current understanding of galac- 
tic chemical evolution or the stellar evolution of low mass stars as they 
pertain to 3 He. 
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